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Boundedness in a chemotaxis-haptotaxis model 

with nonlinear diffusion 
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This article deals with an initial-boundary value problem for the coupled chemotaxis- 
haptotaxis system with nonlinear diffusion 

ut = S/ ■ {D{u)'Vu) — xV • (uVu) — ■ (uS/w) + /ru(l — u — w), x G ft, t > 0, 

< vt = Av — V + u, X G n, t > 0, 

wt = —vw X G n, t > 0, 

K 

under homogeneous Neumann boundary conditions in a bounded smooth domain 
0 C M”, n = 2,3,4, where xA A* given nonnegative parameters. The dif- 
fusivity D{u) is assumed to satisfy D{u) > 6u^~^ for all u > 0 with some 5 > 0. 

It is proved that for sufficiently regular initial data global bounded solutions exist 
whenever m > 2 — For the case of non-degenerate diffusion (i.e. D{0) > 0) the 
solutions are classical; for the case of possibly degenerate diffusion (T)(0) > 0), the 
existence of bounded weak solutions is shown. 

Math Subject Classification (2010): 35K55 (primary), 35B40, 92C17 
Keywords: global existence, boundedness, nonlinear diffusion, chemotaxis, hapto- 
taxis 


1 Introduction 

Already in early stages of cancer, malignant tumours may possess the ability to invade tissue in 
their neighbourhood, with harmful effects on the organism ([23]). Among the many mathemat¬ 
ical models that have been developed for the description of the progress of cancer in different 
stages (see, for instance, [4, 9, 25, 12, 32, 8] and the refercences therein), in [9] Chaplain and 
Lolas introduced the following chemotaxis-haptotaxis system as a model describing the process 
of cancer invasion; 


ut = DAu — V • {xuVv) — V ■ {^uVw) + — u — w), 

X € Cl, 

t > 0, 


TVt = Av — V + u, 

X € Cl, 

t > 0, 

(1.1) 

Wt = —vw + r]w{l — u — w), 

X € Cl, 

t > 0. 
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Herein, L), x, A* > 0, ry > 0 and r G {0,1} are given parameters and u denotes the density 
of cancer cells, v represents the concentration of so-called matrix degrading enzymes (MDEs) 
and w is used to describe the non-diffusive concentration of the extracellular matrix (ECM). 
In addition to their undirected random movement, cancer cells are attracted by MDE and 
macromolecules that are bound in the ECM. Accordingly, their motion is biased toward higher 
concentrations of these substances with this process being known as chemotaxis, if the attractant 
is diffusible, and as haptotaxis for attractants fixed in the tissue. In addition, cancerous cells 
reproduce and compete amongst themselves and with healthy cells for space and nutrients, so 
that source terms of logistic type arise in the first equation. Matrix degrading enzymes, which 
are produced by cancer cells, diffuse, decay, and decompose the tissue around them for more 
living space. A choice of ry > 0 embodies the ability of the ECM to remodel back to a normal 
level. However, since degradation occurs much faster than the re-establishment of the tissue, 
a choice of ry = 0 seems justified, and, indeed, is taken in many of the mathematical works 
concerning systems like (I.I), including the present one. 

If all effects of the extracellular matrix are ignored by letting w = 0 and if additionally /j = 0, 
(1.1) turns into the famous Keller-Segel model [18] for chemotaxis, which has brought forth a 
large amount of mathematical literature during the last decades. We refer to the survey articles 
[15, 14, 3]. In this context one point of particular interest is the question whether solutions 
exist that blow-up in finite time or whether all solutions are global and, maybe, even bounded. 
And in fact, blow-up is known to occur and, according to more recent results, to be a generic 
phenomenon if n > 2 or if n = 2 and the initial mass uq is sufficiently large ([17, 13, 57, 27]). 
The presence of logistic terms (yx > 0, but still re = 0) not only leads to colourful dynamics (cf. 
[19, 31, 30]), but presumably supports the global solvability as well. Nevertheless, proofs of this 
belief, both for the parabolic-elliptic (r = 0) variant [49] and for the parabolic-parabolic (r = 1) 
counterpart [55], require the restriction to the two-dimensional setting or that of sufficiently large 
values of yi > 0. Although there are some closely related models where solutions are global (cf. 
[2, 5, 35, 43, 59, 58]) and although global weak solutions are known to exist ([21]), which, under 
certain conditions on the parameters, become smooth after some time ([21]), it is still unknown 
whether finite-time blow-up may occur if n > 3 and /i > 0 is small. 

Returning to models featuring a nontrivial third component, let us first report on some results 
on haptotaxis-only models, i.e. models with x = 0. For x = /^ = ?? = 0, in [28], the local 
existence and uniqueness of classical solutions have been shown. In [26], a similar model (with 
nonlinear kinetics of the ECM, in that +u in the second equation was replaced by +uw [32]) 
with X = 0 = ^ but yi > 0 was proven to have global bounded weak solutions; the existence 
of a global classical solution was established in [48]. Asymptotic properties of solutions were 
investigated in [24]. For another related model (see [1]), which additionally involves an equation 
governing the evolution of oxygen or other nutrients, in [50] global existence was shown. 

Also for a haptotaxis system including the effect of remodeling tissue (ry > 0) global classical 
solutions are known to exist, at least under the condition that the logistic growth is strong if 
compared to the effects of haptotaxis and regrowth, more precisely, if yr > ^ry [36, 10]. 

As to models involving both chemotactic and haptotactic effects, there have been results guaran¬ 
teeing global existence of classical solutions for ry = 0 under the condition that the chemotactic 
and haptotactic densities depend on the value of u in such a way that x{'^)u aiid ^{u)u are 
bounded [38, 10]. 

For system (1.1) with instantaneous diffusion of enzymes and without tissue remodeling (i.e. 
r = 0 = ry), Tao and Wang [40] proved global existence of solutions for n = 2 or n = 3 and 
sufficiently large yr. Tao and Winkler [45] gave an explicit condition on yr to ensure global 
existence and boundedness and investigated attraction of solutions to the steady state (1,1,0). 
In [46], they improved the condition on yr, so that it coincides with the best one known for 
the parabolic-elliptic Keller-Segel system with logistic source ([49]), and also gave an explicit 
smallness condition on wq, under which w asymptotically becomes negligible. Existence and 
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uniqueness of a global classical solution to ( 1 . 1 ) with 77 > 0 and r = 0 have been obtained in 
[47] in the spatially two-dimensional setting. 

The existence of global solutions to the fully parabolic system (1.1) (that is, r = 1) has been 
established in [39] under the condition of 77 = 0 and 77 being sufficiently large as compared to y. 
Recently, Tao [37] and Cao [ 6 ] proved the boundedness of solutions in two- or three-dimensional 
domains, respectively. 

The models mentioned above described the random part of the motion of cancer cells by linear 
diffusion. As pointed out in [34, Section 6 ], more appropriately from a physical point-of-view, 
however, it might be considered as movement in a porous medium. We will therefore assume D 
to be a nonlinear function of u and shall deal with the following chemotaxis-haptotaxis system 


ut =V • {D{u)Vu) — yV • (uVv) — ■ (uVw) -j- / 2 u(l — u — w), 

vt =Av — V + u, 
wt = — vw, 
du dw dv 
du dv dv ’ 

u{x,0) =uo{x), v{x,0) = vo{x), 777 ( 2 ;, 0) = 7770 ( 3 ;) 


X G D, 

t > 

0, 

(1.2) 

X G D, 

t > 

0, 

(1.3) 

X G D, 

t > 

0, 

(1.4) 

X G cID, 

t > 

0, 

(1.5) 


X G 

n 

(1.6) 


in a bounded domain C M”' with smooth boundary, n G {2,3,4}. The functions uq, vq, wq are 
supposed to satisfy the smoothness assumptions 

77o € C''^(i2) with Tto > 0 in and uq ^ 0, 

< 77o G iy^’°°(fl) with 77o > 0 in Q, ^ 

7770 G for some a > 0 with 7779 > 0 in H and = 0 on dQ. 


Due to the condition = 0 it is already ensured that ^ = 0 on dO, and for t > 0, so that 
(1.5) is equivalent to 


du 
dv 

We furthermore assume that 


dv ^ dw dv 
XUt, -= — = 0 


dv 


dv dv 


X G dQ, t > 0. 


D G (^^([O, 00 )) and D{u) > 6u^ ^ for all 77 > 0 


( 1 . 8 ) 


with some 5 > 0 and tti > 1 which will be subject to additional assumptions. In addition, if 
D{u) satisfies 


D(0) > 0, 


(1.9) 


the diffusion is non-degenerate and one may hope for classical solutions. 

For certain values of m, in [39] global existence of classical solutions to (1.2)-(1.6) was established 
under these conditions. Nonetheless, the question of boundedness of solutions to system (1.2)- 
(1.6) is still open. Thus it is meaningful to analyze the following question: 

Which size of m is sufficient to ensure boundedness of solutions to (1.2)-(1.6)? (Ql) 

Variants of the quasilinear Keller-Segel system with logistic source have been studied in [22, 5, 
51, 52, 7, 60] for instance. Inter alia, Li and Xiang [22] showed global existence of bounded 
solutions to a fully parabolic quasilinear Keller-Segel system with logistic source and degenerate 
diffusion, that is, without assumption (1.9). Hence another natural question to ask is; 

Do bounded solutions to (1.2)-(1.6) with degenerate diffusion exist? (Q2) 

It is our goal in this work to give answers to (Ql) and (Q2). Namely, for non-degenerate and 
degenerate diffusion both, we will show the existence of global-in-time solutions to (1.2)-(1.6) 
that are uniformly bounded. Our main results read as follows: 
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Theorem 1.1. Let n G {2,3,4} and Ll C M"" be a bounded domain with smooth boundary and 
X,i,h > 0- Suppose that the initial data {uo,vo,wo) satisfy (1.7). Then for any 5 > 0 and 

m>2-- (1.10) 

n 

there is C > such that for any function D satisfying (1.8) and (1.9), system (1.2)-(1.6) has a 
classical solution {u,v,w) € (C''^(li x [0,oo) D x (0,oo))))^ which exists globally in time 

and satisfies 

for all t G (0, oo). 

In the case of possibly degenerate diffusion, that is without the assumption (1.9), (1.2) — (1.6) 
possesses at least one global bounded weak solution: 

Theorem 1.2. Let n G {2,3,4} and 11 C M”' be a bounded domain with smooth boundary and 
X,^,h > 0. Suppose that the initial data {uo,vo,wo) satisfy (1.7). Then for any <5 > 0 and 
m> 2 — ^ for any function D satisfying (1.8), system (1.2)-(1.6) has a weak solution {u,v,w) 
in the sense of Definition f.l below that exists globally in time and is bounded in the sense that 
(1.11) holds. 

Remark 1.3. In the case n = 2, condition (1.10) reads m > 1 and coincides with the condition 
posed in order to obtain the global existence result in [41]. For n = 3, (1.10) is stronger than 
the condition m > H in [41] and we leave open the question here whether the solutions to 
(1.2)-(1.6) for Iy < m < I are bounded. 

Remark 1.4. The chemotaxis-haptotaxis system therefore has bounded solutions under the 
same condition on m as the pure chemotaxis system with w = 0 without logistic source [42]. 
For fi = 0 this condition is essentially optimal, cf. [56]. 

The rest of this paper is organized as follows. In the following section, we will prepare the 
later proofs by collecting some useful estimates. Section 3 is devoted to deriving a differential 
inequality for the quantity + Jq | Vu]^'^, p, > 1, whose boundedness is asserted in Lemma 
3.9 and which can be used to prove the boundedness of solutions and thus Theorem 1.1. Relying 
on the existence of classical solutions in the non-degenerate case, in Section 4 we will then 
complete the proof of Theorem 1.2 by an approximation procedure. 

2 Preliminaries 

In this section we collect a few short, but helpful results, most of which have already been proven 
elsewhere. In some cases, for the convenience of the reader, we take the liberty of stating the 
conditions differently from the original source, in such a way that the lemmata are still covered 
by the original proof, but become more easily applicable in the present situation. For example. 
Lemma 2.7 below was proven in [37] for a system of PDFs differing from (1.2)-(1.6) by the first 
equation. The only property of u required in the proof, however, is nonnegativity, so that this 
difference plays no role, and we state the lemma accordingly. 

Let us begin with an estimate for a particular boundary integral that enables us to cover possibly 
non-convex domains. 

Lemma 2.1. Let 12 C M"' be a bounded domain with smooth boundary, let q G [l,oo) and 
M > 0. Then for any rj > 0 there is Crj > 0 such that for any v G (7^(12) satisfying ^ = 0 on 
dLl and /q |Vu| < M the inequality 

f <r? / |V|Vun2 + C^ 

holds. 
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Proof. This has been proven in the course of the proof of [16, Prop. 3.2] by a combination 
of an estimate of on the boundary in terms of the curvature of dQ and jVup with the 

embedding ^ a fractional Gagliardo-Nirenberg inequality 

and Young’s inequality (cf. Lemmata 2.2 - 2.5 and estimate (3.10) of [16]). □ 


We will make use of the Poincare inequality in the following form: 


Lemma 2.2. Let fl C M” he a bounded smooth domain. Let a > 0. Then there exists C > 0 
such that 


\u 


W1.2(Q) ^ C 


llV^r]] 


L2(n) 


+ 



for all u G IT^’^(fl). 


Proof. Assuming to the contrary that there exists a sequence (un)neN C ly^’^(H) such that 
ll^n|lL2(f^) = 1 and llun||i2(Q) > n(]lVunll^2(Q) + (/f^bnr)“) for any n G N, we can find a 
subsequence thereof converging to some v G weakly in IT^’^(n) and strongly in L^(n) 

so that ]|u||j;^2(f2) = 1 . Due to ju|“ < lim inf„_,.ool^n|“ < hminf„_,.oo(i)“ = 0, at the same 
time we obtain u = 0, a contradiction. □ 


Also the Gagliardo-Nirenberg inequality will be used in a less common version (for a similar 
variant see [53, Lemma 3.2]): 

Lemma 2.3. Let D C M” be a bounded smooth domain. Let r>l,0<q<p< oo, s > 0 be 
such that 

1 1 1 

( 2 . 1 ) 


1 1 1 

— < —I—. 
r n p 


Then there exists c > 0 such that 


\u\\LP{n) < c 


L’’(C) 


I^IIl9(o) + ll^^l 


L=(C) 


for alluG n L'J{n), 


where 


a = 


1 _l_ i _ i 
q r 


1 • 


( 2 . 2 ) 


Proof. Compared to the standard version of the Gagliardo-Nirenberg inequality ([29, p.l26], [11, 
Thm. 10.1]), we include the possibility of p, g G (0,1), but restrict ourselves to the case oi q < p. 
The conditions q < p and y < ^ ^ ensure that a G [0,1]. For a proof of the case 1 < p < g, we 




refer to the aforementioned theorems of [29, 11]. For g < 1 < p, let c := -j —j-- Then c G [0,1] and 

q p 

^ = 1. Holder’s inequality with exponents ^ and yields llttll 2 ,i(Q) < II^IIlp(q) 

and we conclude from the usual version of the Gagliardo-Nirenberg inequality that 


I 

\L^{n) 


I II II 11^ II 11^—^ II 11^ II 11^^^ — 

m\LP(n) - '-"1 iFllwi.’-(n) ItIIwi.’-(q) ItIIlp(q) 


(l-c)(l-b) 


I'“IIl9(o) 


l-i 

for b = and with some Ci > 0 so that 


(l-c)(l-6) 


II ^ ^ II II 1—c(i—fe) II II 1—c(i—6) ^ II \\a II ||1 —a 

0Ilp(C 2) — *-^2 IPIIwi.’'(n) ~ ^2 iFllrvi.r(n) II'W'IIl9(q) 


where C 2 = and 

b 


1 - 


a = 


I — c + cb 


1 + i _ 1 1-1 1-1 

^ _P I 2. 


1 1 

q p 


1-^ 
P 


TZX ■ i+±-i 

q p nr 


l + i_ 1 + 1 _1 

n r q 
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coincides with the expression for a given in (2.2). For q < p < 1 let c := j G [0,1]. Then 


1 

p 


cp + = 1 and by Holder’s inequality with exponents cp and we have 


|c(l—6)+l—c 


'“'IIlp(Q) — ll^llLi(n) ll“llL'2(f2) — ^3 ll'w|lvKi'’'(n) ll'“llL'?(n) 


= c3\\u\ru 


VKi’’-(r 2 ) 


1^11 


1—a 

Li(n) 




with C 3 from the usual Gagliardo-Nirenberg inequality, where b = —r and thus a := cb = 

q~^ n r 

i_i 

1 ■ In both cases, in ||n||^p(m < max {( 72 , Cs} ||n||^i,rm\ we may employ the 

q'n r 

Poincare inequality (Lemma 2.2) to obtain Cp > 0, so that with (74 = max {(72, ( 73 } • Cp we 
achieve < <74 ||Vu||2.(q) + (74 ||n||}^7(n) ||u||2.(f^), where ||u|| 2 .(f 2 ) < 

C 5 ||n|| 2 ^s(Q) if s > q with some C 5 obtained from Holder’s inequality. If, on the other hand 

s < q, then l|n||/,s(Q) < Cq < Cq ll'*^llLs(f 2 ) — 2 II^IIlp(C 2 ) ^7 

with Cq > 0 and C-/ > 0 from Holder’s and Young’s inequality, respectively, and d = {- — 


-)/(- — -). 
qJ/ \s p' 


□ 


When we have to estimate powers of norms, we will often withont notice combine the Gagliardo- 
Nirenberg inequality with the following elementary estimate. 

Lemma 2.4. For every a > 0 there is (7 > 0 such that for all x,y a (0, 00 ) we have 

(x + yT < G(a:“+ ?/“). 


Proof, {x + ?/)“ < (2 max {x, y})" < 2" max {x“, y"} < 2"(x“ + y“). 


□ 


Next let us give a basic property of the total mass of cancer cells that can be checked easily. 


Lemma 2.5. LetTyO, letu € C‘^’^{Qx{0,T))r\C^{Qx[0,T)) solve (1.2) with some functions 
v,w G (7^7(f2 X (0, T) n(7°(II X [0, T)), D € (7^(M) and 0- Uu w are nonnegative, 

u satisfies 



X < m* := max 



(2.3) 


for all t € ( 0 , T). 


Proof. Thanks to the homogeneous boundary condition, we directly integrate (1.2) with respect 
to space. Using the nonnegativity of u and w and the Gauchy-Schwarz inequality, we obtain 


d 

dt 


I u = fj, —yL / W 

In Jn Jn 




/ uw u — y 

Jo. Jo. Jo, 


< h 




2 


on (0,T). With an ODE comparison argument, this leads to (2.3). 


□ 


The following properties, tnrning information on the first solution component into boundedness 
assertions about the second, can be derived by invoking the variation-of-constants formula for 
V and estimates for the heat semigroup. 

Lemma 2.6. Let T G (0,oo], let v G C^’^(Q x (0, T)) n C^(fi x [0,T)) solve (1.3) with some 
function 0 <u G C^{Ll x [0,T)) and let M > 0. 

(i) Assume that ||w(-T)||ii(f 2 ) < M for all t G (0,T). Then for any s G [1,;^), there exists 
c = c{M, s) > 0 such that 


v{-,t)\\L‘>{n) < c for all t G (0,r). 
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(ii) Assume that ||w(-, < M for all t G (0,T). Then for any s G [1, there exists 

c = c{M, s) > 0 such that 


l|V?;(-,t)||L»(n) < c for all t G (0,T). 


(in) Assume that p > | and ||'w(-, < M for allt G (0,T). Then there is c = c{p,M) > 0 

such that 

ll^(•>^)llL°°( 0 ) ^ c for all t G (0,r). 

(iv) Assume that p > n and < M for all t G (0,T). Then there is c = c{p, M) > 0 

such that 

||V?;(-,t)||^oo(s^) < c for all t G (0,r). 

Proof, (i) Let s G [1,;^). We denote by the (Neumann-) heat semigroup and use its 
positivity to estimate 




L^{n) — 




L°°( 0 ) 


L'>( 0 ) 


LA^) Jo 


i,“(r 2 ) 


-h [ +u{-,t)) 

Jo 

+ Jo \\u{-,T)\\Ls^n^dT 

^)llLi(o) 

ll^o|lLoc(f,) + J^ (ci(l + {t- dr 


for all t G (0, T), where ci is the constant from the L^-L^-estimate in [54, Lemma 1.3 (i)], Ai > 0 
is the first positive eigenvalue of —A in fl, and u{t) = f^u(x,T)dx, r G (0, T), denotes the 
spatial average of u. Thanks to s G [1,;^), we have —f(l — -j) > —1 and hence the last 
expression is bounded independently of t G (0, T). 

(ii) Let s G [1, From [54, Lemma 1.3 (iii) and Lemma 1.3 (ii)] we obtain C 2 > 0 and 

C 3 > 0, respectively, such that 


l|Vu(-,t)|Us(Q) < C2||Ve*(^ + [ l|Ve(* ^^u{-,T)\\Ls^Q)dT 

Jo 

nf 

< C2||Vuo||Loc(f^)+ C3 / (l + (t-r)"^"?(^"-))e“'^'(*"^)||u(-,r)||i,i(Q)(ir 

Jo 

for any t G (0, T). Again, the inequality —^ — ^(1 — i) > —1 and the boundedness assumption 
on u provide us with a time-independent bound for this expression. 

(iii) Let p > ^. Lemma [54, Lemma 1.3 (i)] gives ci > 0 such that 


\vit)\\ 


L°°{n) — 




^ IbollLoo(Q) + ^ ^ci(l-I-(t 

< ll^o|lLoo(f2) + ^ (^Ci(l-|-(t 


+ / 


Jo 



-r) 2p) ||u(.,t)- n(-,T)||^p(f^)e + e ^M|M(r 

-r)-i)2Me-^^(^-") +e-(*-")|L!|-i |k(-, t)||^i(^)) dr 



dr 


for every t G (0, T), which again results in boundedness due to —^ > —1. 

(iv) Let p > n. Passing to the limit p —)• oo in [54, Lemma 1.3 (iii)] and from a variant of 
[54, Lemma 1.3 (iv)] (given in the form needed here e.g. in [20, Lemma 3.1]), similarly to the 
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previous part of the proof, we obtain C 2 > 0 , C 4 > 0 such that 


llVull 


L°°(Q) — 


Ve*(^-b 


^^0 


^ Jo 


Ve(*-r)('^-i)u(-,r) 


dr 




<C 2 ||Vuo||ioo(t^) + C4y^ (l + (t-r) 2 2p)e 




\LP{n) 


dr 


for any t G (0,T), which again, because —^ ~ ^ > “Ij entails uniform-in-time boundedness of 

Due to the fact that (1.4) is an ODE, w can be represented explicitly in terms of v. Following 
an observation from [37], this provides a one-sided pointwise estimate for —Aw: 


Lemma 2.7. Assume that x > 0, ^ > 0 and /x > 0, a > 0, T > 0, and let the nonnegative 
functions v,w G C'^’^(D x (0, T)) n 0*^(0 x [0, T)) solve (1.3) and (1.4), respectively, with some 
nonnegative u G C^(Q x [0,T)) and let w{-,0) G C'^"''“(D) be positive in D. Then 

—Aw{x, t) < ||rc(', 0)||/,oo(Q) • v{x, t) + K for all x gTL and t G (0, T), 


where 

K := ||Au;(-,0)||ioo(o) +4||vV<;W||i.e(f,) + (2.4) 

Proof. This can be found in [37, Lemma 2.2]. □ 

Let us finally recall the following local existence resnlt, which has been established in [41] by 
means of a standard fixed point argument. 

Lemma 2.8. (Local existence) Let n G {2,3,4} and D C M” be a bounded domain with smooth 
boundary. Let > 0 assume that D satisfies (1.8)-(1.9), and the initial data fulfils 

(1.7). Then there exists a maximal time Tm a^ G (0, 00 ] and a triple {u,v,w) of functions from 
C^{0, X [0, Tmax) nC'^’^(D X {0,Tmax))) solving (1.2)-(1.6) classically inLlx (0,Tmax) o.nd such 
that 


either Tmax = 00 or limsup (||u(-, t)||ioo(f 2 )-h llu(-, t)||tyi,cx>(Q)-h ||r(;(-, t)||i 4 /i,oo(Q)) = oo.(2.5) 

t/^Tmax 

Moreover, 


u > 0, u > 0, and 0 <w < |ltt)ollL°°(n) in Q x (0,Tmax). (2.6) 

Proof. See [41, Lemma 2.1], □ 

3 Boundedness of solutions. Proof of Theorem 1.1 

The purpose of this section is to prove Theorem 1.1, which is concerned with the non-degenerate 
case. 

For the rest of this section let us fix n € {2,3,4} and the initial data satisfying (1.7). Given 
X) ? > 0 and a fnnction D, by (u, v, w) we will always denote the corresponding solution to 
(1.2)-(1.6) given by Lemma 2.8 and by T^ax its maximal time of existence. The value of K will 
be as defined by (2.4). 

Now we proceed to establish the main step towards our boundedness proof. Motivated by [42], 
we establish a differential inequality for the expression u^P |Vu|^'^ on ( 0 , T^ax) for p,q > 1 , 
from which we will be able to derive a bonnd on P /n 
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Lemma 3.1. Let > 0, 5 > 0 and m> 1. For any function D satisfying (1.8) and (1.9) 

and for any p > 1 the solution {u,v,w) to (1.2)-(1.6) satisfies 


1 d 


.,P ^ 


+ 


pdt 

■ X {P ~ ^) f m+l 


— I D{u)uP-^\Vu\^ + 


(3.1) 


<- 


[ uP-™+'|V^P + elKllL-=(f,) [ uPv + {p + ^K) [ uP-p 

JQ. JQ. Jo, 


,p+i 


for all t G (OjTjnax) and with K as in (2.4). 

Proof. Multiplying (1.2) by and integrating over 12 we obtain 

f u^ + {p—l) f D{u)u^~^\Vu\‘^ < x(p “ 1) [ • Vu + ^(p — 1) f u^~^S/u-S/w 

P dt Jq Jq Jq 


+p uP-p uP+^ -p uPw 

Jo. Jq, Jq. 


(3.2) 


on (OjTmax). Here we split the integral containing D into one part we will keep until we need 
it in Section 4 and one part that can be used to cancel less favorable contributions by other 
integrals: According to (1.8), 

(p-1) [ D{u)uP-^\Vu\^ [ D{u)uP-^\Vuf+ [ u^+P-^\Vu\^, (3.3) 

Jn ^ Jn 2 Jq 

which again holds on the whole time-interval {0,Tjnax)- In the next integral in (3.2) an appli¬ 
cation of Young’s inequality yields 

X(p-l) [ uP-^Vu-Vv < [ u^+P-^\Vuf + ~ [ uP-^+^\Vv\'^ (3.4) 

Jn ^ Jn 0 Jn 

on (0,Tmaa;). According to an integration by parts and Lemma 2.7, we also have 

?(P- 1 ) 


?(P-1) [ uP ^ 
Jn 


Vtt • Vw = 


< 


P 

fjp- 1 ) 

P 


f uP{—Aw) 

Jn 

[ uP{\\wo\\L^in)-v + K) 

Jn 


< ?II«^o||l“(o) f uPy + fK f uP on(0,r„ 
Jn J n 

with K as in (2.4). Substituting (3.4), (3.3) and (3.5) into (3.2), we obtain 


(3.5) 


1 d 
p dt 

< 


[ uP + ^—- [ L»(u)uP-2|Vup + f u 

Jn Jn 4 Jn 

[ uP-^+^\\/v\‘^ + C\\wo\\L^^n) [ uPy 

0 Jn Jn 

+ ^K [ uP + p [ yP - p [ uP+^ -p [ yPw 

Jn J Jn Jn 


P+™-3|Vup 


on (0, Tjyihx 


)■ 


Finally inserting uP^^ 3|Vy|2 — ^ |2 using the nonnegativity of f^yPw, 

we arrive at (3.1). □ 

One of the terms on the right-hand side of (3.1) that has to be dealt with is the integral uPy. 
Since at the moment high powers of u and y are out of reach for our estimates, we will use the 
Gagliardo-Nirenberg inequality and gradient terms to control this integral: 
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Lemma 3.2. Let > 0- Let 5 > 0 and m > 2 — For all p G [l,C!o) and r] > 0 there 

exists a constant C = C{p,r]) > 0 such that for any function D fulfilling (1.8) and (1.9) the 
solution to ( 1 . 2 )-( 1 . 6 ) satisfies 


p+m —1 


Vu 2 (.ji) 


L2(Q) 


Clko||L->(Q) / uP{-,t)v{-,t) <r] 

Jn 

for all t € ( 0 ,rmax). 

Proof. Because (n — 4)p < 0 < 2(m — 1) and m > 2 — L guarantee that 
p + m — 1 n — 2 


+ C 


(3.6) 


2 pf 


> 


2 n 


and 


P- I 


it is possible to fix 7 € ( 0 , 1 ) so small that still 

p + m — 1 n — 2 




> 


2-7 


2 n 


and 


p+m —1 

n —2 

2 

2 n 

p - 

2-7 

n 

p+m —1 

n —2 

2 

2 n 


< 2 , 


< 2. 


(3.7) 


We set r = and r = and note that by Lemma 2.6(i) applied to s = there exists 

C > 0 such that ||f (•, t)|Lr+o^ - ^ ^ ^ {0,Tjnax)- An application of Holder’s inequality 

'1 j' ^j 

then asserts 

[ u^{-,t)v{-,f) < ( [ u^^(-,t)^ ■( [ u’’(-,t)V for alH G (OjTmax). 

Jn \Jn / \Jn J 

As 0 < ^q:A_ < and, by (3.7), \ < ^ + from the Gagliardo-Nirenberg inequality 

and Lemma 2.5, we obtain constants ci,C 2 > 0 such that 


pr 


p+m— 1 

= \\u 2 


2p 


< 


< 


I p+m —1 
I 2pr 

LP+m-\ (f2) 

p + m-1 -- ^ , ^P 


/ ,, P + m —1 ,,—^ 

Cl llVa—III+) 


^—^-a,, p+m— 1 ,,— p.^.m—1 —j-^ 

1-1 I — 2 - ||P+m-l' -I I p+m-1 


+ C 2 


_ P+m—1 

Vtt 2 


2p 


P+m— 1 

L 2 ( 0 ) 


+ 1 


2 + ll'W ^ II 2 

J^p+m —1 ^ p+m —1 


on ( 0 , Tpfiax ) 


with 


o = 


p+m— 1 p+m— 1 1 2—7 

2 2 pr p m 1 pn 


2 pr 

p+m —1 I 1 1 

2 n 2 


p+m —1 _ n —2 ' 
2 2 n 


Since 


2p 


p _j_ — 1 ^ p+m —1 _ n —2 


p 


2-7 


< 2 , 


by (3.7), an application of Young’s inequality thus completes the proof: Given 7 > 0 it provides 
C > 0 such that 




\\'Wo\\L--in) J^uP{-,t)v{-,t) 


<p 


p+m-1 

Vu 2 (.^t) 


p+m-1 , 

Vu 2 

2p 

p+m —1 

+2(0) 

2 

+ c 

L2(0) 



+ 1 


for all t G (0,Tmaa:). 


□ 
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The next part of the differential inequality we are searching for is given by the following lemma. 


Lemma 3.3. Let ^ 0- Let 5 > 0 and m > 1. For q € [l,oo) there exists a constant 

C = C{q) > 0 such that for any function D with (1.8) and (1.9), the solution to (1.2)-(1.6) 
obeys 


IfL 

q dt 






+ 


[ iVlVul^p 
Jn 



(3.8) 


on (0, TjYidx ) • 

Proof. This can be proven by merging the estimate from [16, Prop. 3.2] with some of the 
calculations from the proof of [42, Lemma 3.3]: We first observe that 

[ \Vv\‘^‘^ = 2 [ jVuj^^-^VuVAu - 2 [ • Vu + 2 [ \Vv\'^‘^-^Vv ■ Vu (3.9) 

Q dt Jq Jq Jq Jq 

holds throughout ( 0 ,Tmaa:)) and use the pointwise identities 

2VvVAv = A\Vv\^-2\D^v\^ and V\Vv\^^-^ = {q - l)\Vv\‘^'^-'^V\Vv\^ (3.10) 

together with an integration by parts to obtain 


i2 |2 

^ v\ 


2 [ iVul^-J-^VuVAu = [ -2 [ 

Jn Jn Jn 

= -{q-l) [ |Vu| 2 ''-^|VlVul 2 | 2 + [ 

Jn Jdn 

-2 [ \Vv\^^-^\D\\^ on{0,Tmax). 

Jn 

Another integration by parts, (3.10) and Young’s inequality make it possible to estimate the 
rightmost term in (3.9) according to 

2 [ \Vv\‘^'i-^VvVu = -2{q-l) [ ulVul^'J-^VjVuj^ • Vu - 2 [ ujVuj^^-^Au 

Jn Jn Jn 

^~-f \Vvf'i-^\V\Vvf\'^ + 2{q-l) [ u^\Vvf 

Jn J n 


<- 


\2q-2 


+ - / ^ on( 0 ,r^,,,), 

n Jn ^ Jn 

where we have used jAu] in the second-last integral. Adding these estimates, we 

arrive at 


“W f |Vu|2^ + 2 [ |Vup''<-^-^ [ |Vu|2^"^lV|Vu| 
Q dt Jn Jn 2 Jq 


2|2 


+ 


( 2 (, - 1 ) + f 

2 Jn Jdn ou 


on (0,Tmax) and thus at (3.8) if we take into account that j 

and use Lemma 2.1, which is applicable by a combination of Lemma 2.6 (ii) and Lemma 2.5, 
to gain c > 0 such that 


[ /“ jVlVun^ + c. on( 0 ,T^„,) 

Jdn Q Jn 


□ 
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The terms on the right-hand side of (3.1) we have not yet treated are +{^K+fi) vP—^ 

For positive /r they could easily be estimated from above by a constant. In order to also cover 
the case of /i = 0 , we prepare the following estimate instead: 

Lemma 3.4. Let > 0. Let S > 0 and m > 1 — Then for any p > 1 and rj > 0 there 

is a constant C = C{r]) > 0 such that for any function D with (1.8) and (1.9), the solution to 
( 1 . 2 )-( 1 . 6 ) satisfies 


f < 7 ] f 
Ju Jq 


\S/u 




+ C 


on ( 0 , TjYiax ) • 


Proof. Once more, an application of Lemma 2.5 shows that f^u(-,t) < ci for some ci > 0. 


2p 


> —r and i < - 

p+m—1 2 n 


Obviously, ^ 2 - n ' 2 ' 2 p 

Nirenberg inequality asserts the existence of C 2 > 0 such that 

2 v 


I ~ n therefore the Gagliardo- 


II 

p+m—1 

U 2 

In 



Lp+m-1 (Q) 


<C 2 


/ 

P+m—1 

,a 

P+m—1 

( 

Vtt 2 

L2(0) 

U 2 


2p 


’ ( 1 -a) 

P+m—1 ' ' 

2 -|- 
Lp+m-1 (Q) 


P+m—1 

U 2 


2p 


P+m—1 
Lp +^-1 ( q ) 


. (1—a) 

P+m—1 ^ •' 


<C2C[ 

on { 0 ,Tmax), where 


[ |Vu 

Jn 


p+m—1 ,o 

2 r 


p+m—1 


2p 


+ C 2 Cf 


P+m—1 


a = 


p+m—1 
2 


p+m —1 

2 p 


p+m—1 I 1 1 

2 w n 2 


and 


P 


-a = 


P 


{p + m-l){l - i) 


p — 1 


p-l-m —1 p-l-m —1 p-l-m— 1 - 1 -^ —1 p—l-|-(m— 1 -|- + 
so that an application of Young’s inequality gives the desired conclusion. 


< 1 


□ 


What we have achieved with the previous estimates is the following: 

Corollary 3.5. Let p,> 0. Let h > 0 and m > 2 — For any 1 < p,q < oo there 
exists a constant C = C{p,q) > 0 such that for any D fulfilling (1.8) and (1.9), the solution to 
( 1 . 2 )-( 1 . 6 ) satisfies 


dt 


I [ uP+ [ [ D{u)uP-^ 

{Jn Jn J P JQ 

+ 2 q [ + f |V|Vu|5p<C f 

Jn Q Jn Jn 


|Vu|2 + 

^p-m+l 


6 p{p - 1 ) 


, p+m — 1 

|Vu 2 


2 {p + m- 1)2 

|Vup + C [ + C (3.11) 

Jn 


on the whole time-interval ( 0 . Tma^ ). 


Proof We employ Lemma 3.2 and Lemma 3.4 with p = or rj = , 

respectively, as well as Lemma 3.1 and Lemma 3.3 and add the resulting inequalities. □ 


In a similar manner as in [42], we deal with the two terms on the right hand side of (3.11). We 
will have a closer look at the conditions on p, q in Lemma 3.8 afterwards. 

Lemma 3.6. Let p, y, ^ > 0, h > 0 and m > 2 — Let p,q > 1 be such that p > m — and 


1 > 


q p — m-\- 1 
q — Ip + m — 1 


^(p -|- m 


_ nq-n+2 
_ g(p-m+l) 

(p -|- m — l)n -|- 2 



(3.12) 


12 























































Then for any r/ > 0, one can find a constant C = C{T],p,q) > 0 such that for any D fulfilling 
(1.8) and (1.9), the solution to (1.2)-(1.6) satisfies 


[ < rj [ + / |V|V?;|«p + (:7 

J 17 J Q. J 17 


(3.13) 


for all t G ( 0 ,rmax)- 


Proof. Because p > we have nqp > mnq — n + 2 and thus nq — n + 2 < nq+npq — nqm = 


- T " . Furthermore, > 1 so that 

p+m—l — (p+m—1) (nq—n+2) ’ ^ — 


< 


nq{p — m + 1), which results in 0 < 
n — 2 — 2 q < 0 and hence — 2 pq-\-p[n — 2 ) < 0 < ng(m— 1 ) + (g —l)(n — 2 )(m — 1 ), which results in 
g(n — 2 )(p —m+ 1 ) = pqn — 2 pq — q{n — 2 ){m — \) <pqn—p{n — 2 )+nq[m—\) — {n — 2 )[m—l) = 
(ng - n + 2 )(p + m - 1 ) and thus \-\ = ^< • 

All in all, by these estimates and (3.12), it is possible to choose 7 > 0 in such a way that 


p — m + 1 
p + m — 1 


(p + m — 1 ) 


_ ng-n+ 2-7 N 
q(p-m+lj 

{p + m — l)n + 2 — n 


< 


g- 1 


(3.14) 


and 


1 1 

-< 

2 n 


{nq — n + 2 — 7 )(p + m — 1 ) 
2 nq{p — m + 1 ) 


and ng — n + 2 — 7 < ng(p — m + 1). (3.15) 


Our main purpose in introducing 7 is to avoid n — 2 in the denominator in the 2-dimensional 
case. For n G {3,4} also using 7 = 0 would be possible. Invoking Holder’s inequality with the 
exponents n- 2 +'y 


[ < ( 

J 17 V-j 


nq—n+2—7 
ng(p —m+1) \ nq 

'll nq — n+ 2 —y 


n 


n — 2+7' 
2 nq nq 

I V'Uj ^-2+7 


on (0, Trnax ) ■ 


1 o .^n 

From the Sobolev embedding VF^’^(n) L"- 2 + 7 (ll) and Lemma 2.2 we obtain ci > 0 and 

C 2 > 0 , respectively, such that for arbitrary sq € [Ij we have 


2nq 


I Vul >^-2+7 


n — 2+7 
nq 


=iiivunr^ <ciiiivurii^,,, 

Ln-2+j (fl) 


<C2 ( ||V|Vuni2,(^) + |||Vun|% j < C2 ||V|Vur||2,(^) + C3 


on (OjTmax), where is the contant given by Lemma 2.6 (ii) corresponding to sq. 

Thanks to the first inequality in (3.15), the Gagliardo-Nirenberg inequality therupon provides 
us with C 4 > 0 such that 


/ 

Jn 


nq(p —m+1) 
'll nq-n+2 —7 


nq—n+2—7 
nq 


2(p-m+l) 

H p+m—1 —■. - .+■ 

11^2 - P+m -1 

^ I I 2nq(p —m+1) 


{p+m-l)(nq-n+2 —i') 


2a(p —m+1) 


< c4(||v«“^ii7+7- 


2(1 —a,)(p —m+1) 
I p+m —1 

I 2 

Lp+m-l (Q) 


+ lk 


2(p-m+l) 

P+m —1 , , —I- ^ 

^ ^ - p+m —1 

^ I I 2 

Lp+m-l (Q) 


(3.16) 


on { 0 ,Tmax), where 


a = 


p+m—l _ (p+m—1)(nq—n+2—-y) 
2 2nq(p—m+l) 

p+m—1 I _ 1 
2 w n 2 
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Thanks to the boundedness of u by Lemma 2.5, we may continue (3.16) to estimate 


nq —n+2 —'Y 
Tiq(p —m+1) \ nq 

'll nq — n+ 2 —f 


p+m —1 


2 a(p—m+1) 
p+m —1 


<C 4 ||Vm 2 +C5 On( 0,rmaa;) 


with appropriate C 5 > 0 obtained from Lemma 2.5. 

Given r? > 0, Young’s inequality gives ce > 0 such that, in conclusion. 


[ < C2C4 ( [ |V|Vr;|''p 

Ju \Jq 


^ p —m+1 
p+m -1 r, \ p+m -1 


'n 


+ C 3 C 4 / |Vrt 2 


p+m -1 \ P+m -1 


p —m+1 

+ C2C5 ( / iVlVul'^p )" +C3C5 

Jq 

p —m+1 

p+m -1 


<C 2 C 4 i^j IVlVul-^P 

+ ? / |V|Vu |«|2 + ^ [ iVn'^p + ce on( 0 ,T^a+. 
Jn Jn 

Furthermore, by (3.14), we also have 


(3.17) 


a{p — m+ 1 ) (? — 1 

p + m — 1 q 

Hence in the first term on the right hand side of (3.17) we can employ Young’s inequality to 
obtain C 7 > 0 such that 


C2C4( / IVlVul^l^ 


|Vm 


p —m+1 

p+m-1 riX^P+Ti-l V f .rri'i V f i, , ^+"1-1,0 

<- / |V|Vr;|''r + 7( / |Vu^~r + C7 


to finally transform (3.17) into (3.13). 


□ 


We can treat the second term on the right-hand side of (3.11) similarly as the first one in Lemma 
3.6: 

Lemma 3.7. Let p, x, ^ > 0, 5 > 0 and m > 2 — Let q > 1 and p > max{l, 2 q^^n -2 “ ^ + 1} 
be such that 

p+m—1 (p+m—l)(2g+n—2) 

2 4nq 


1 2 
q p + m — 1 


p+m—1 I 1 1 

2 n 2 


(3.18) 


Then for any r/ > 0 there exists a constant C > 0 such that for any D with (1.8) and (1.9) the 
solution to ( 1 . 2 )-( 1 . 6 ) fulfils 


[ < T] [ + [ iVlVul^p + G 


(3.19) 


for all t € ( 0 ,Tmax)- 


Proof. Since (n — 2 ){q — 1) > 0 implies 2 q+n -2 — (3-18), we can find 7 > 0 such 


that 


nq 

2q+n-2-7(q-l) 


> 1 and 


p + m — 1 


p+m—1 _ (p+m—l)(2g+n—2—7(g—l)) 
2 4nq 

p+m—1 I 1 1 

2 n 2 


< 


(3.20) 
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Using Holder’s inequality with exponents 
yields 

[ < 

Jn 


__ 

2q+n-2-'y{q-l) 


and 


nq 


nq 


nq—2q—n+2+^{q—l) (n— 2 + 7 )((j—1) 


_ 2 ng_ 

2q+n-2- j(q-l) 


2q + n-2 - r(q-l) (q-l)(n-2+j) 

nq ! r 2nq \ nq 

I |V'U|”“^+T' 


on ( 0 , Trr 


where we may again use Sobolev’s and Poincare’s inequalities to obtain 


(<;-l)(n- 2 + 7 ) 

2nq \ nq 2q^ 

|Vu|^^ =|||Vu|''||| 

/ 2,n-2+7 (Q) 


< ciiiivuniX 


in) 

2q-2 


< C 2 ||V|Vur|| L^{n) + |||Vu| 


2 

<?|| Q 


£a 
L 9 


for all t G {0,Tmax) with sq G [1, We have 2nq > nq + nq > 2q + n > 2q + n — 2 


so 


that 0 < 


< 


4nq 


n ' 4nq 

C 3 > 0 with 


— p+m—l — (p+m—l)(2g+n—2) ' 

i)( 2 g 

Anq 


Moreover, - m + 1 < p entails i = i + %^ = 


2q^n^~2'^ /o | o\ 

- + and hence the Gagliardo-Nirenberg inequality yields 


2nq 


||2g+n-2-7(g-l) 


2 g+n-2—fCg —1) 
nq 


,, p+m—l ,, —;- 7 

= +2 - 


4ng 


I ||„ p+m-l 

< C3 l|Vtt 2 


(p+m-l){2q+n-2 —/(ij-l)) (f7) 
p+m 

U 2 


4b , , 4(1-6) 

I-r p+m—l ,,—) -+ 

p+m—l - n - p+m — l 


IL 2 (Q) 


I, p+m-l 

+ ||U 2 


I P+m—l 

I 2 

L p+m-l (O) 


P+m — l 


46 


< c 4 (|| Vu “||£"+1 


L p+m-l (n) 


on (0, T„ 


where we have used Lemma 2.5 to find C 4 > 0 and where 

b = 


p+m—l _ (p+m—l)(2g+n—2— 7 (g—l)) 
2 Anq 


p+m—l I 1 1 

2 w n 2 


Here (3.20) ensures that p^^_i < | and concluding as in the previous proof, by twofold 

application of Young’s inequality (with exponents in the first step) we obtain 

(3.19). □ 


The conditions in the previous lemmata involve some assumptions on p and q that are not 
immediately seen to be simultaneously satishable. With this lemma we ensure that they are. 
In its proof we rely on the fact that m > 2 — + 

Lemma 3.8. Let m > 2 — There exist unbounded sequences {pk)keii^ iQk)keii such that for 
each A: G N 


qk >1 
Pk > max 


and Pk > 1 
2qk{n - 2) 


1 >- 


2 qk + n — 2 
qk Pk-m+l 


qk-lpk + m-l 
2 


— m + 1, m — 


{pk + m- 1) 


n — 2 


nqk 


_ nqk-n+2 n 

gkiVk-m+T) 

{pk + m — l)n + 2 — n 


1 

qk Pk + m-1 


Pfc+m-l _ (pfc+m-l)(2qfc+n-2) 
2 _ Anqk 

Pfc+m-1 _|_ 1 _ 1 


and hence pk and qk fulfil all requirements of p and q from Lemmata 3.6 and 3.7. 


(3.21) 

(3.22) 

(3.23) 

(3.24) 
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Proof. In a first step we will show that it is possible to choose pk, Qk in such a way that the 
conditions 


f 2 qk{n-2) n-l n-2 

Pk >max<^ --^-- -m + l,2qk - m + l,m - 

[2qk + n-2 n nqk 

Pk <{‘^qk - l)(ni - 1) + — 

n 


(3.25) 

(3.26) 


are satished, which evidently already imply (3.22), but, in a second step, can also be used to 
derive (3.23) and (3.24) as well. 

There is Qo > 0 such that for every q > Qo 


max 


2?(ra - 2) , ^ 2 q{n - 1) n-2 

— m + 1,- m+ l,m - 


2 q + n — 2 


n 


nq 


n — 1 

= 2q - m + 1 


n 


as can be seen by considering the limit of each of the expressions as g —>• oo. By the condition 
on m we know that m + ^ — 2 > 0 and hence it is possible to hnd Qi > 0 such that for every 
q > Qi we have 

2 q\m-\ -2) >m—1 


n 


and hence 


2 q 2q , 2q 2q(n — 1) 

2 qm — m — 2q + 1 ->-h 2 q, i.e. {2q — l)(m — 1) H-> -, 

n n n n 

so that it becomes possible to choose unbounded sequences fulhlling (3.21), (3.25) and (3.26) 
simultaneously. This is already sufficient for them to satisfy (3.23) and (3.24) also. Indeed, 
(3.26) ensures that 


npk < {2qk — l)(m — l)n + 2 qk = 2qk + 2 nmqk — 2nqk — mn + n 


and hence 


npqk + nmqk — 2nqk + 2 qk — npk — mn + 2n — 2 
>npkqk + nmqk — “^^nqk + 2 qk — {2qk + 2 nmqk — 2nqk — mn + n) — mn + 2n — 2 
>npkqk - nmqk + n-2, 


so that 
1 > 


npkqk - nmqk + n-2 


_ _ nqkjpk -m+1) - {nqk - n + 2) 

npkqk + nmqk - 2nqk + 2 qk - npk - mn + 2n - 2 {qk - l)(?^^'fc + nm - 2n + 2) 

n(7fc—rx+2 

/ V — - 

{Pk + m- 1)- 


_ qk Pk-ni + l 
qk-'^Pk + m-l 

Moreover, (3.25) entails 


qk{Pk-m+l) 


“^qk- 


n — 1 


n 


{pk + m — l)n + 2 — n j 


+ 1 — m < Pk, i.e. 2 qk{n — 1) + n — mn < npk 


and thus 


2 nqk — 2qk — n + 2 < npk + mn — 2n + 2, 


which shows that 

2 nqk -2 qk-n + 2 


1 > 


_ 2qt,+n-2 
qk 2n 


‘^qk 


npk + mn — n + 2 — n Pk+m i j_ _ i 


+ i i Pk + m-1 


Pfc+m-l _ {pk+m-l)(2qk+n-2) 
2 __ 

Pk+m-1 I _ 1 
2 n 2 


□ 
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Having made sure that there are (arbitrarily large) exponents, for which the lemmata before 
can be applied, we now use them in order to obtain the following boundedness result. 

Lemma 3.9. Let x,^ > 0, 5 > 0 and m > 2 — Let p,q G (l,oo). Then one can find a 
constant C > 0 depending onp,q such that for any function D satisfying (1.8) and (1.9), we 
have 

[ uP{;t)+ [ |Vu|25(-,t) <C7 (3.27) 

J Q J Q. 


for all t G (0,Tmax)- 

Proof. Observing L^{Ll) ^ L^{Ll) for r > s, we may without loss of generality enlarge p and q 
and do so in a manner that all conditions on p and q listed in Lemma 3.6 and 3.7 are satisfied, 
which is possible due to Lemma 3.8. An application of Corollary 3.5 provides us with constants 
Cl, C2, C3, C4 > 0 such that 

[ uP+ [ |Vu|25^+ci + [ iVup^ + cs [ \V\Vv\'^f 

dt \Jq, Jn J Jn Jn Jn 

< C4 f + C4 f + C4 

Jo. Jq 


holds on {0,Tjnax)- Invoking Lemma 3.6 and Lemma 3.7 upon the choice of p 
we find C5 > 0 such that 




fL 

dt 



(3.28) 


on {0,Tmax)- Moreover, as a consequence of Lemma 3.4 there is cg > 0 fulfilling 


C2 


[ uP f \Vu 

Jn ^ Jn 


2 \ +CQ 


(3.29) 


and combining (3.28) and (3.29) shows that 




< C 5 + C 6 


for all t G {0,Tmax)- If we define y{t) := f^uP(-,t) + |Vi;(-,t)p'i', t G (0,Tmax), and C7 := 

C5 + C6, the inequality reads 


y'(t) + C 2 y(t) < C 7 for all t G (0, T^ax)- 


Upon an ODE comparison, this entails boundedness of y and hence (3.27). □ 

A direct consequence is the following assertion on boundedness: 

Corollary 3.10. Let p, x, ^ > 0, 6 > 0 and m > 2 — For any p G [1, 00) there is C > 0 such 
that for any function D obeying (1.8) and (1-9), the solution to (1.2)-(1.6) fulfills 

\\u{-,t)\\LP{n) <C and \\Vv{-,t)\\^p < C 

for any t G {0,Tmax)- 

Since p = 00 is not covered by this corollary, we set out to improve the norms which can be 
controlled. In the case of v nothing more than a short application of —L'^-estimates is needed: 
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Lemma 3.11. Let x, > 0, 5 > 0 and m > 2 — There is C > 0 such that for any function 
D obeying (1.8) and (1.9), the solution to (1.2)-(1.6) satisfies 


lb(■)^ C* 

for all t € [ 0 ,Tmax)- 

Proof. We employ Corollary 3.10 for some p > nto gain a time-uniform bound on ||u(-,f)||^p(Q^. 
Then Lemma 2.6 (hi), (iv) show that there is ci > 0 such that ||Vu(-, < ci and 

ll^(•)*)llL°°( 0 ) ^ Cl for all t G {0,Tmax)- □ 

Turning Corollary 3.10 into boundedness of u is more difficult. The lack of a bound on Vu; makes 
||n(t)ll^oo(Q) inaccessible for general boundedness results like [42, Lemma A.l]. We proceed 
somewhat similarly to [37]: 

Lemma 3.12. Let > 0, h > 0 and m > 2 — There exists C > 0 such that for any 

function D satisfying (1.8) and (1.9), the solution to (1.2)-(1.6) fulfils 


l'*^(') II ^ 


for all t G [ 0 ,Tjnax)- 

Proof. From Lemma 3.11 we obtain ci > 0 such that ||Vu(-, t)||^oo(Q) < ci and ||t'(',i)||icx>(fi) < 
Cl for all t G {0,Tmax)- 

First let us fix pQ > |(m — 1) so large that G |) for all p G [po; oo)- 

From Lemma 3.1 we see that 


1 d 

pdtj^ (p + m-l)'^ 




£±!f^|2 ^ G(P-I) /■ „,-„+! 


f uP-^+^\Vv\ 

0 Jn 

+ ?lho|lLoc(m / uPv + {p + f,K) I uP - P I 
Jn Jn Jn 


w 


p+i 


with K as in (2.4) and for all p > po,t G {0,Tmax)- We may conclude 






|Vn^^2 P <p^C2 I w 


<p C2 / W 

In 


[ + C3P [ uP-pp [ uP+^ (3.30) 

Jn Jn Jn 

[ + C3P [ uP - pp [ uP~^^ +p^C2|n| 

Jq Jq. 


for all p > po and all t G (0, Tmax), where 


2 X 


C2 = cf— and C3 = f, \\wo\\/ q) ci + p + f,K + 1 and 5 <- 


6 Ap{p — 1) 


§ II “-U II 1^.20 (11) I 

According to Young’s inequality we may estimate 


2 {p + m — ly 


C3p [ uP <pp [ uP+^ + cf^^lJ^I- 

Jn Jn P +1 V P / 

Furthermore, with a = have ^ = (5 “ ^)o + so that the Gagliardo-Nirenberg 

inequality gives 04,05 >0 such that by Young’s inequality 


h^n) <CA{\\VT\\ltl^) 11^112^(0))' + C5 Il7^llii(ii) 

(5 

L'^(n) 


0 /ki2 


2 p^C2 


l|V<^||? 


, , , 2 2 n p^C 2^5 

+ ( C5 H-rT:(— TF— 2 


n -|- 2 n -|- 2 S 


2 

P(n) 
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for all (/9 G 

p+m—1 

If we apply this to (p = u 2 and insert these estimates into (3.30), we arrive at 


A 

dt 


/“'+/ 


/ < 


C3P 

P+l 


p +1 


p P\n\+p C 2 C 5 + -—( 


2n p^C 2 R 


re + 2 n + 2 6 


)2C4 


p+m— 1 

u 2 


+p^C2|fI| 
(3.31) 

for any p > po and t G (0, Tmax)- 

We then recursively define pk ■= 2pk-i + 1 — m for A; G N and note that po > |(rre — 1) and 
m>lhy an inductive argument lead to po>^ > + 2~^){m — 1) > 0 for all fc G N (at least 

if m > 1; if m = 1, pfc = 2 ^pq), which obviously entails the existence of cg, cy > 0 such that 


^ Pk . 

^6 < ^ < C7 

for any A; G N. Furthermore, we fix T G (0, Tmax) and denote 


(3.32) 


Mk ■= sup 
te(o 


p f uP'^{t). 
,T) Jn 


From (3.31) we hence conclude 

f < f Uq*" 

Jq J q. 

+ j + (^P^C2C5 + ret5-tp^+2^ Mk-i +p^C2|0|^ dr. 


that is 


Mk < ||reo||fp,(f^) +p +PfcC2|ll| + C8Pfc+^M|_4 


for all A: G N with 


, 2 C 2 C 4 . n + 2 n n 

C 8 :=C 2 C 5 + (——) 2 n^d 2 . 
re + 2 


Either there is a subsequence {pki)i such that 

cspt"Ml_, < 


LPk(n)+P +PfcC2|ll| 


and hence 


|'w|lLoo(f2) <limsupM^^' 


1 


Z—>-oo 


<limsup (2 [||reo||^^p^f2 )+/i +p^C2|12 |l)< ||reo|| 


Z—>-oo 


\L°°{n) 


+ 21 + 1 


or for all sufficiently large k we have Mk < 2csp^~^‘^ M^_^ and hence 

Mk < c^pl^^Ml_^ < C9C7(2'^)"+2m|_i 
for some cg > 0 and for all A: G N by (3.32) and thus 

Mk < a^M^_i for all A; G N with a = max {cgcy, 1} 2”'“'"^. 
Induction leads to 


(3.33) 
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Since by another inductive argument we have k + 2-^(^ — j) < for all A: € N and 

because is bounded, we hence obtain from (3.33) that 

_L 2'' — 1 

iV4^" < < {a^MoY-e =: cio- 

Taking fc —>■ 00, we finally arrive at 


sup ||u(-,t)l4ocm) < CIO 
te(o,T) 

and due to the arbitrarity of T, this shows the claim. □ 

We are now in position to pass to the proof of Theorem 1.1. 

Proof of Theorem 1.1. According to Lemma 2.8, there exists a classical local-in-time solution 
{u,v,w) to (1.2)-(1.6) on [0, T^ax)- By Lemmata 3.12, 3.11 and 2.8, u, v and w are uniformly 
bounded on [0,Tmax) and consequently (2.5) makes it possible to conclude T^ax = 00 . □ 


4 Degenerate difFusion and global weak solutions. Proof of 
Theorem 1.2 


The goal of this section is to prove Theorem 1.2. In the absence of (1.9), the first equation 
of system (1.2)-(1.6) may be degenerate at u = 0 and we cannot hope for classical solutions. 
Therefore we introduce the definition of weak solution to (1.2)-(1.6) which we shall pursue here. 

Definition 4.1. Let T G (0, oo] and 12 C ffi"" be a bounded domain with smooth boundary. 
A triple {u, v, w) of nonnegative functions defined on x (0, T) is called a weak solution to 
(1.2)-(1.6) on [0,r) if 

(i) u G LL([0,T);L40)), V G LL([0,T); W^-^II)) and w G LL([0,T); 

(ii) D{u)VuGLlY[0,T);LYn)y, 

(in) {u,v,w) satisfies (1.2)-(1.6) in the sense that for every ip G C'“(n x [0,T)), 


/ uift- uo{x)ip{x,0)d- 
Jo Jn Jn 


X = — 


f [ D{u)Vu-Vip + x [ [ uVv-Vip 

Jo Jn Jo Jn 


+C / / u'^w-Vip + p / / uip — p / / Yip 

Jo Jn Jo Jn Jo Jn 

-fi / uwip] (4.1) 


holds as well as 

rT 


[ [ vipt- [ vo{x)p{x,0) = - [ [ Vv -Vif - [ [ v(p+ [ [ uif (4.2) 

Jo Jn Jn Jo Jn Jo Jn Jo Jn 


and 


In particular, if T 
( 1 . 2 )-( 1 . 6 ). 




wo{x)ip{x,0) 



(4.3) 


oo can be taken, then (rt, v, w) is called a global-in-time weak solution to 
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In order to obtain a weak solution to (1.2)-(1.6), we start with the approximate problem, for 
e G (0,1) given by 


Urf = V • {D^{Ue)VUe) - yV • (UeVVe) - • {UeVWe) + - We), 


Vet = AUe -Ve + Ue, 


< Wet = -VeWe, 

Oue _ dvE _ dwE _ f) 

diy du du ’ 

rt£(x,0) = uq { x ), U£(x,0) = vq { x ), r(;£(x,0) = r(;o(x) 


X G 11, t > 0, 
X G 11, t > 0, 
X G H, t > 0, 
X G 511, 1 > 0 
X G H 


where De{s) = D{s + e) for all s > 0. 

Since D G (^^([0, oo)) fulfils D{s) > for all s > 0, also for De the estimate 

De{s) = D{s + e)> 5{s + (4.5) 


holds with the same values of <5 and m, which, notably, are independent of e. Furthermore, 

De{0) = D{e) > > 0. (4.6) 

Accordingly, Theorem 1.1 and the lemmata from its proof become applicable so as to yield 
global classical bounded solutions {ue,Ve,We) to (4.4) and the following bounds, uniformly in 
e > 0; 

Lemma 4.1. Let n G {2, 3,4}, y, ^ > 0 and suppose that m > 2 — ^ and <5 > 0. Assume that 
the initial data {uq,vq,wq) satisfy (1.7) and the diffusion function D fulfills (1.8). Then for any 
p > 1 there exists a constant C > 0 such that 


ll'^e(') IIl°°(Q) —C*! 

(4.7) 

Il'^e(‘5^) AC*, 

(4.8) 

K'e(') L°°{n) AC*, 

(4.9) 

r [ De{Ue)uP-^\VUef <C{l+t), 

Jo Jn 

(4.10) 

[ [ uT+^-^\VUef <C{l+t), 

Jo Jn 

(4.11) 

[ \VVe{;t)\^ <C, 

Jn 

(4.12) 

[ [ \VWef <C{l+t). 

Jo Jn 

(4.13) 


for all t G (0, oo). 


Proof. The bounds in (4.7), (4.8) and (4.9) are direct conclusions of Theorem 1.1. The estimates 
(4.10), (4.11) can be gained directly from Lemma 3.1 by integration with respect to time. From 
(4.8) we furthermore know there exists a constant fulfilling (4.12). Moreover, according to 
Lemma 2.7, we have 


/ |Vu;£|^ = - / WeAWe < \\woe\\Loo / VeWe + (K + 1) / We, 

J Q, 

so that the boundedness of in (4.13) results from the bounds on Ve and We 


□ 


We can use these bounds in a straightforward manner to extract convergent subsequences of 

Ve,We-. 


(4.4) 
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Lemma 4.2. There exist v,w € oo), VL^’^(r2))nL°°(Ox (0, oo)) and a sequence {£k)ke'N \ 

0 such that 

in lL((0,oo),l2(1^)), (4.14) 

in L°°{Q X {0,oo)) (4-15) 

in L°°{^ X (0,oo)) (4.16) 

in Lf^^{{0,oo),L‘^{n), (4.17) 

in LUiO,oo),L\n), (4.18) 

in L°°(fl X (0, oo)) (4.19) 

as e = Ek \ 0- 

Proof. This directly follows from the estimates (4.8), (4.9), (4.13). □ 

Due to the nonlinearities involved, we require better convergence properties of u^. As prepara¬ 
tion for an Anbin-Lions-argument, let us first ensure boundedness of the time-derivatives in a 
certain weak sense. The reasoning is similar as in [44, Proof of Theorem 1.1] and [52, Proof of 
Theorem 1.3]. 

Lemma 4.3. Suppose that m > 2 — Let 6 > maxjl, y}. Then there is r > 1 such that for 
any T > 0 there exists C > 0 such that 



Vvs Vu 


Ws ^ w 

VWe Vw 
Ws ^ W 


i4)i 


Lp(o,Ty,(wf’^(n)r) 


< C 


(4.20) 


holds true for any e G (0,1). 

Proof. We choose r > 1 so large that Wq’^(D) 1T^’°°(D). Using that (^“(D) is dense in 

IUq’^(D), we may express the norm as 


{<) 









sup 

CeCo-(0),||C||^y2.r-(n)<i 



(4.21) 


In order to estimate this integral, let C & be snch that ||C|Ivk 2 .'-(o) — ^ ^*4 hence 

IjVClIi^oc < Cl and ||Clli oo(Q) < Cl with Cl being the embedding constant of Wq’^(D) ‘- 5- 1U^’°°(D). 
We then multiply the first equation in (4.4) by u^~^C npon an integration by parts, obtain 

I [ {u%-C = -{ 6 - 1 ) [ D,{ue)ut-^\Vue\\- [ D,(u,)uf-'Vu, • VC 
0 Jn Jn Jn 

+X{6 - 1) / U^e~^S/Ue ■ S/VsC + X [ ulVVe ■ VC 
Jo. J Q 

+c(0 - 1) / 4 -^Vue • Vu;,C + C / u^We • VC 
J Q, J D. 

+p [ u^^{l -Ue- We)C =: h + ... + Ij (4.22) 

Jn 

on {O^Tmax)- By (4.7) and (4.9), there are 02,03 > 0 such that ||w£|li;,oo(Qx(oT)) — ^2 and 
II^£|Il°°(Ox(o,t)) — ^3’ that on ( 0 -,TfYiax) 

Ihl < tk J^ul{l-Ue-We)\C\ < AiC2-(l-FC2-bC3)-||CllL°°(Q) 1^1 ^ /iC2(l-bC2-hC3)ci |D| =: C4. (4.23) 
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Also /4 can be estimated rather easily. With the aid of (4.12), we find C5 > 0 such that 
|Vn£(-,t)p < C5 for all t € (0,T) and hence 

U 4 I < XC 2 f I^VeVCI < XC 2 f f IVVeP / |VCp) <XC 2 c|ci=:C 6 On(0,T^aa:)- (4.24) 
J \J S7 J Q. / 

We do not have as convenient bounds for Iq, but we can prepare an estimate of its integral 
by means of (4.13). In order to do so, we use Young’s inequality to see that with c^ ;= 
^C2Ci max {1, |n|} 


1^61= ? [ • VCI < ?C2Ci [ (iVtCep + 1) < C7 + C7 [ on(0,Tmax)- (4.25) 

Jq Jq. Jq 

We want to handle the terms containing |Vu£p by means of (4.11). Therefore, let us fix p > 1 
such that 9 > max |p, | ^ Then by Young’s inequality 


\h\=xio-i) 


j 4-‘ 

Jn 


VUe ■ VVsC 




<X(9 - 1) ^ |Vu,|2) . IICIUco(o) 

<x(9 - 1) + c^) ci = cg / + cg (4.26) 

\ Jn J Jn 


on {0,Tmax), where we have set cg = x(0 — l)c2^ ^ cg = x(0 — l)c5Ci. Again by Young’s 

inequality, we obtain 


l^5|=^(0-l) 


^S/Ue ■ VWeC 


<C(0-1) ( / uf-2|Vue|"+ / iVn;, 


'n 


\L<=°{n) 


<?(« -1) (ii”=iif-(L7„T)) I ^ I VU..I") ■ iiciIl»(o), 

<.e(0-l) [ u^+P-^lVuep + f |Vr 

V Jn Jn 


'We\ Cl 


<C 10 [ nr+P-3|Vn,|2 + Cii [ \VWef on {0,Tmax), 
Jn Jn 


(4.27) 


with the abbreviations cio = ^{9 — and cn = ^{9 — l)ci. For the remaining two 

terms, we note that D^{u^) = D{u£ + e) < maxo<s<c2+i D{s) =: d and that, by our choice of 9, 
29 — p — m + 1 > 0, so that 


\h\ = 


De{Ue)ul ^VUe • VC 


< 


< D,{u,)vP-^\VUe\^ + Cl 

on {0,Tmax) and 


( [ D,iu,)uP-‘^\Vus\^ + [ Il,(n,)nf-^|Vn,|2 V IIVCIIl-(o) 
\Jn Jn J 

(4.28) 


|/i| = (0-1) 




<{9- l)||u,||J^Cf(^) (^Il,(u,)nr2|Vu, 


< C12 / De{Ue)ul ‘^\VUe 


(4.29) 
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for any t G (0, Tmax), if we set cu = (^ — l)c2 ^ci. Finally combining (4.23) - (4.29) with (4.21) 
and (4.22), we obtain 




< 


2 —I (-fl + • • • + 

Li((0,T);(Vy2-’'(O))*) Jo 

Tr rT, 

29—p—m+l 
2 


(Cl2+Cl) [ [ Ds{Us)uP ‘^\S/Ue\‘^ + cidcf ^ 

JoJn JoJn 


+ (cg + Cio) / / ^IVUgp + (cii + Cy) / / + (cg + Cg + Cy + C4)T. 


0 


0 


Here we can apply (4.10), (4.11) twice (also for p = 2) and (4.13), so that we obtain an e- 
independent bound for the right-hand side and thus have proven (4.20). □ 


Lemma 4.4. Let m > 2 — Then there exists a subsequence (ej)jgN of the sequence provided 
by Lemma f.2 such that for all T > 0 


Ue ^ U 

a.e. in D X (0, T), 

(4.30) 

Ue ^ U 

in L2((0,r),L2(F!)), 

(4.31) 

D{u)Vu 

in L\i0,T),L\n)). 

(4.32) 


Proof. We fix p > 1 and 9 = Given T > 0, the bound from (4.11) then shows that there 

is C > 0 such that 1^ < C* for any e G (0,1) and together with the uniform bound from 

(4.7) we infer that (uf)e6(o,i) is bounded in L^(0, T; PF^’^(n)). From Lemma 4.3, we know that 
((Ue)t)ee(o,i) is bounded in Li((0,T); (IFo^’^(H))*). Since W^’^(n) L^(n) ^ (Wo’^(fl))*, 

we therefore conclude from the Aubin-Lions compactness lemma [33, Cor. 4] that (rtf)£6(o,i) 
is a relatively compact subset of the space L^(0,T; L^(Q)), whence there exists a sequence of 
numbers s = Sj '\ 0 such that —)■ in L^(0, T; L^(fl)) along this sequence with some 
function in this space, which, along a subsequence, ensures (4.30). Together with (4.7), 
Lebesgue’s dominated convergence theorem therefore asserts (4.31) as well. 

The boundedness of as warranted by (4.7) entails the uniform L°°-boundedness of so 

that inserting p = 2 into (4.10) already ensures boundedness of i7e(tte)Vue in T^((0, T), L^(H)), 
and after extracting a further subsequence we may assume that D^{ue)^Ue converges weakly 
in L^((0, T), L^(n)). In the interest of identification of the limit, we resort to a primitive of 
say Gefx) := D^{.s)ds = D(s)ds. We already have observed that V[G£(tie)] = D^{u^)S/u^ 
is weakly convergent in L^((0, T), L^(H)). Moreover, by (4.30) + e ^ u a.e. in H x (0,T) so 

that 

rUe{-)+e rui-) 

G,{ue{-))= / Dis)ds^ / D{s)ds =: G{u{-)) 

Jo Jo 

a.e. in Ll x (0,T) and thus, again as consequence of (4.7) and Lebesgue’s dominated convergence 
theorem, in L^((0, T), L^(fl)). Therefore the weak limit of V[Ge(u£)] has to coincide with 
VG(u) = D{u)Vu, so that we arrive at (4.32). □ 


Proof of Theorem 1.2. The convergence properties asserted in Lemma 4.4 and Lemma 4.2 are 
sufficient to pass to the limit in each of the integrals making up the corresponding weak for¬ 
mulation associated with (4.4) so that {u,v,w) is a weak solution to (1.2)-(1.6) in the sense of 
Definition 4.1. The boundedness of {u,v,w) in the sense of (1.11) results from (4.30), (4.15), 
(4.19) and (4.16). □ 


5 Acknowledgements 

The hrst author has been supported by China Scholarship Council (No.201406090072) and in 
part by National Natural Science Foundation of China (No.11171063). 


24 







References 


[1] A. R. Anderson. A hybrid mathematical model of solid tumour invasion; the importance 
of cell adhesion. Math. Med. Biol, 22(2):163-186, 2005. 

[2] K. Baghaei and M. Hesaaraki. Global existence and boundedness of classical solutions for a 
chemotaxis model with logistic source. C. R. Math. Acad. Sci. Paris, 351(15 - 16):585-591, 
2013. 

[3] N. Bellomo, A. Bellouquid, Y. Tao, and M. Winkler. Toward a mathematical theory of 
keller-segel models of pattern formation in biological tissues. Math. Models Methods Appl. 
Sci., 25(09):1663-1763, 2015. 

[4] N. Bellomo, N. Li, and P. K. Maini. On the foundations of cancer modelling: selected 
topics, speculations, and perspectives. Math. Models Methods Appl. Sci., 18(04):593-646, 
2008. 

[5] X. Cao. Boundedness in a quasilinear parabolic-parabolic Keller-Segel system with logistic 
source. J. Math. Anal. Appl., 412(1):181-188, 2014. 

[6] X. Cao. Boundedness in a three-dimensional chemotaxis-haptotaxis model. arXiv preprint 
arXiv:1501.05383, 2015. 

[7] X. Cao and S. Zheng. Boundedness of solutions to a quasilinear parabolic-elliptic keller- 
segel system with logistic source. Math. Meth. Appl. Sci., 37(15):2326-2330, 2014. 

[8] M. A. Chaplain and A. R. Anderson. Mathematical modelling of tissue invasion. Cancer 
modelling and simulation, pages 269-297, 2003. 

[9] M. A. Chaplain and G. Lolas. Mathematical modelling of cancer invasion of tissue: dynamic 
heterogeneity. NHM, l(3):399-439, 2006. 

[10] J. Fan and K. Zhao. A note on a 3D haptotaxis model of cancer invasion. Appl. Math. 
Res. Express. AMRX, (l):74-86, 2014. 

[11] A. Friedman. Partial differential equations. Holt, Rinehart and Winston, Inc., New York- 
Montreal, Que.-London, 1969. 

[12] R. A. Gatenby and E. T. Gawlinski. A reaction-diffusion model of cancer invasion. Cancer 
research, 56(24):5745-5753, 1996. 

[13] M. A. Herrero and J. J. L. Velazquez. A blow-up mechanism for a chemotaxis model. Ann. 
Scuola Norm. Sup. Pisa Cl. Sci. (4), 24(4):633-683 (1998), 1997. 

[14] T. Hillen and K. J. Painter. A user’s guide to PDE models for chemotaxis. J. Math. Biol., 
58(1-2):183-217, 2009. 

[15] D. Horstmann. From 1970 until present: the Keller-Segel model in chemotaxis and its 
consequences. 1. Jahresber. Deutsch. Math.-Verein., 105(3):103-165, 2003. 

[16] S. Ishida, K. Seki, and T. Yokota. Boundedness in quasilinear keller-segel systems of 
parabolic-parabolic type on non-convex bounded domains. J. Differential Equations, 
256(8);2993-3010, 2014. 

[17] W. Jager and S. Luckhaus. On explosions of solutions to a system of partial differential 
equations modelling chemotaxis. Trans. Amer. Math. Soc., 329(2) :pp. 819-824, 1992. 


25 


[18] E. F. Keller and L. A. Segel. Initiation of slime mold aggregation viewed as an instability. 
J. Theor. Biol, 26(3):399-415, 1970. 

[19] K. Kuto, K. Osaki, T. Sakurai, and T. Tsujikawa. Spatial pattern formation in a 
chemotaxis-diffusion-growth model. Phys. D, 241(19):1629 - 1639, 2012. 

[20] J. Lankeit. Chemotaxis can prevent thresholds on population density. Discrete Contin. 
Dyn. Syst. Ser. B, 20(5):1499-1527, 2015. 

[21] J. Lankeit. Eventual smoothness and asymptotics in a three-dimensional chemotaxis system 
with logistic source. J. Differential Equations, 258(4):1158-1191, 2015. 

[22] X. Li and Z. Xiang. Boundedness in quasilinear Keller-Segel equations with nonlinear 
sensitivity and logistic source. Discrete Contin. Dyn. Syst., 35(8):3503-3531, 2015. 

[23] L. A. Liotta and T. Clair. Cancer: checkpoint for invasion. Nature, 405(6784):287-288, 

2000 . 

[24] G. Li^canu and C. Morales-Rodrigo. Asymptotic behavior of global solutions to a model 
of cell invasion. Math. Models Methods Appl. Sci., 20(09):1721~1758, 2010. 

[25] J. S. Lowengrub, H. B. Frieboes, F. Jin, Y.-L. Chuang, X. Li, P. Macklin, S. M. Wise, and 
V. Cristini. Nonlinear modelling of cancer: bridging the gap between cells and tumours. 
Nonlinearity, 23(1):R1-R91, 2010. 

[26] A. Marciniak-Czochra and M. Ptashnyk. Boundedness of solutions of a haptotaxis model. 
Math. Models Methods Appl. Sci., 20(3):449-476, 2010. 

[27] N. Mizoguchi and M. Winkler. Blow-up in the two-dimensional parabolic Keller-Segel 
system. 2013. Preprint. 

[28] C. Morales-Rodrigo. Local existence and uniqueness of regular solutions in a model of 
tissue invasion by solid tumours. Math. Comput. Modelling, 47(5-6):604-613, 2008. 

[29] L. Nirenberg. On elliptic partial differential equations. Ann. Scuola Norm. Sup. Pisa (3), 
13:115-162, 1959. 

[30] K. Osaki, T. Tsujikawa, A. Yagi, and M. Mimura. Exponential attractor for a chemotaxis- 
growth system of equations. Nonlinear Anal., 51(1):119-144, 2002. 

[31] K. J. Painter and T. Hillen. Spatio-temporal chaos in a chemotaxis model. Phys. D, 
240(4-5):363-375, 2011. 

[32] A. Perumpanani and H. Byrne. Extracellular matrix concentration exerts selection pressure 
on invasive cells. European J. Cancer, 35(8): 1274-1280, 1999. 

[33] J. Simon. Compact sets in the space L^(0,T;R). Ann. Mat. Pura Appl. (4), 146:65-96, 
1987. 

[34] Z. Szymahska, C. Morales, M. Lachowicz, and M. A. J. Chaplain. Mathematical modelling 
of cancer invasion of tissue: the role and effect of nonlocal interactions. Math. Models 
Methods Appl. Sei., 19(2):257-281, 2009. 

[35] Y. Tao. Boundedness in a chemotaxis model with oxygen consumption by bacteria. J. 
Math. Anal. Appl, 381(2):521-529, 2011. 

[36] Y. Tao. Global existence for a haptotaxis model of cancer invasion with tissue remodeling. 
Nonlinear Anal. Real World Appl, 12(l):418-435, 2011. 


26 


[37] Y. Tao. Boundedness in a two-dimensional chemotaxis-haptotaxis system. arXiv preprint 
arXiv:1407.7382, 2014. 

[38] Y. Tao and C. Cui. A density-dependent chemotaxis-haptotaxis system modeling cancer 
invasion. J. Math. Anal. AppL, 367(2) :612-624, 2010. 

[39] Y. Tao and M. Wang. Global solution for a chemotactic-haptotactic model of cancer 
invasion. Nonlinearity, 21(10):2221, 2008. 

[40] Y. Tao and M. Wang. A combined chemotaxis-haptotaxis system: the role of logistic 
source. SIAM J. Math. Anal, 41(4):1533-1558, 2009. 

[41] Y. Tao and M. Winkler. A chemotaxis-haptotaxis model: the roles of nonlinear diffusion 
and logistic source. SIAM J. Math. Anal, 43(2):685-704, 2011. 

[42] Y. Tao and M. Winkler. Boundedness in a quasilinear parabolic-parabolic Keller-Segel 
system with subcritical sensitivity. J. Differential Equations, 252(1):692-715, 2012. 

[43] Y. Tao and M. Winkler. Eventual smoothness and stabilization of large-data solutions in a 
three-dimensional chemotaxis system with consumption of chemoattractant. J. Differential 
Equations, 252(3):2520-2543, 2012. 

[44] Y. Tao and M. Winkler. Locally bounded global solutions in a three-dimensional 
chemotaxis-Stokes system with nonlinear diffusion. Ann. Inst. H. Poincare Anal. Non 
Lineaire, 30(1):157-178, 2013. 

[45] Y. Tao and M. Winkler. Boundedness and stabilization in a multi-dimensional chemotaxis- 
haptotaxis model. Proc. Roy. Soc. Edinburgh Sect. A, 144:1067-1084, 10 2014. 

[46] Y. Tao and M. Winkler. Dominance of chemotaxis in a chemotaxis-haptotaxis model. 
Nonlinearity, 27(6): 1225, 2014. 

[47] Y. Tao and M. Winkler. Energy-type estimates and global solvability in a two-dimensional 
chemotaxis-haptotaxis model with remodeling of non-diffusible attractant. J. Differential 
Equations, 257(3):784-815, 2014. 

[48] Y. Tao and G. Zhu. Global solution to a model of tumor invasion. Appl. Math. Sci. (Ruse), 
l(45-48):2385-2398, 2007. 

[49] J. I. Tello and M. Winkler. A chemotaxis system with logistic source. Comm. Partial 
Differential Equations, 32(6):849-877, 2007. 

[50] G. Walker and G. E. Webb. Global existence of classical solutions for a haptotaxis model. 
SIAM J. Math. Anal, 38(5):1694-1713, 2007. 

[51] L. Wang, Y. Li, and G. Mu. Boundedness in a parabolic-parabolic quasilinear chemotaxis 
system with logistic source. Discrete Contin. Dyn. Syst., 34(2):789-802, 2014. 

[52] L. Wang, G. Mu, and P. Zheng. On a quasilinear parabolic-elliptic chemotaxis system with 
logistic source. J. Differential Equations, 256(5): 1847-1872, 2014. 

[53] M. Winkler. A critical exponent in a degenerate parabolic equation. Math. Methods Appl. 
Sci., 25(ll):911-925, 2002. 

[54] M. Winkler. Aggregation vs. global diffusive behavior in the higher-dimensional Keller- 
Segel model. J. Differential Equations, 248(12):2889-2905, 2010. 


27 


[55] M. Winkler. Boundedness in the higher-dimensional parabolic-parabolic chemotaxis system 
with logistic source. Comm. Partial Differential Equations, 35(8): 1516-1537, 2010. 

[56] M. Winkler. Does a ‘volume-filling effect’ always prevent chemotactic collapse? Math. 
Methods Appl. Sci., 33(l):12-24, 2010. 

[57] M. Winkler. Finite-time blow-up in the higher-dimensional parabolic-parabolic Keller-Segel 
system. J. Math. Pares Appl. (9), 100(5):748-767, 2013. 

[58] T. Xiang. Boundedness and global existence in the higher-dimensional parabolic-parabolic 
chemotaxis system with/without growth source. J. Differential Equations, 258(12):4275- 
4323, 2015. 

[59] C. Yang, X. Cao, Z. Jiang, and S. Zheng. Boundedness in a quasilinear fully parabolic 
Keller-Segel system of higher dimension with logistic source. ArXiv e-prints. Mar. 2015. 

[60] J. Zheng. Boundedness of solutions to a quasilinear parabolic-elliptic keller-segel system 
with logistic source. J. Differential Equations, 259(1):120 - 140, 2015. 


28 


